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We present the hrst application of statistical spectroscopy to study the root- 
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■ mean-square value of the parity nonconserving (PNC) interaction matrix el- 
ement M determined experimentally by scattering longitudinally polarized 
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■ neutrons from compound nuclei. Our effective PNC interaction consists of a 
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standard two-body meson-exchange piece and a "doorway" term to account 
for 0~ spin- flip excitations. Strength functions are calculated using realis- 
tic single-particle energies and a residual strong interaction adjusted to fit 



^ . the experimental density of states for the targets, 238 U for A ~ 230 and 

104,105, 106.108PJ £ or a ~ 100. Using the standard Desplanques, Donoghue, 
and Holstein estimates of the weak PNC meson-nucleon coupling constants, 
we find that M is about a factor of 3 smaller than the experimental value 
for 238 U and about a factor of 1.7 smaller for Pd. The significance of this 
result for refining the empirical determination of the weak coupling constants 
is discussed. 
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I. INTRODUCTION 



Studies of the weak interaction in strongly interacting systems have been of great interest 
in nuclear physics for a number of reasons. First of all, on the experimental side, observables 
may be obtained cleanly from measurements because: (1) the symmetries respected by 
the weak interaction are different from those respected by the strong interaction, allowing 
the weak interaction sector to be uniquely disentangled from the purely strong interaction 
sector; and (2) the strength of the weak interaction is sufficiently small that its matrix 
element between definite eigenstates of the strong interaction Hamiltonian may be extracted 
perturbatively. On the theoretical side, knowledge of the weak interaction observables is a 
rich source of insights into the behavior of many-body systems. In those cases where the weak 
interaction is known, the experimental observable provides an opportunity to learn about 
the nuclear many-body wave function, giving an often unique glimpse at a specific feature of 
the eigenstate of the nuclear Hamiltonian. If, on the other hand, the weak interaction is not 
completely known (as is the case for the hadronic weak interaction), the measurements offer 
the possibility of learning about the weak interaction itself. To realize this goal, the main 
theoretical task is to calculate the properties of the nuclear many-body system as accurately 
as possible with a Hamiltonian that is known to describe the initial and final states in the 
matrix element. 

Experimental information likely to shed light on the nature of the hadronic weak in- 
teraction is available in both light and heavy nuclei. Measurements of weak-interaction 
observables in these cases are driving the development of theoretical formulations of the 
many-body problem along two principal lines distinguished by the significance placed on the 
interpretability of individual matrix elements. The shell model formulation of the many- 
body problem is an example of a nonstatistical approach, and here one envisions individual 
matrix elements, such as those of operators describing transitions between discrete nuclear 
levels, to be calculable and hence meaningful. In light nuclei |l| the data have been available 
for many years, and these data are being analyzed using such a shell-model approach. 



2 



Our particular interest is in heavy nuclei. In this case 0, measurements based on lon- 
gitudinally polarized neutron scattering and focusing on the hadronic parity nonconserving 
(PNC) interaction have become available only in the last 5-10 years, and development of 
the microscopic statistical approach is being stimulated by the availability of this data. The 
main focus of the latest experiments || has been the strength function of the effective PNC 
interaction in nuclei in the mass A ~ 230 and mass A ~ 100 regions of the periodic table. 
The root-mean square (RMS) matrix element M of this interaction is the physically mean- 
ingful quantity, as evidenced by the fact that the nuclear spacing between the resonances 
involved is on the order of 10 - 100 eV; note that it is not possible to obtain individual 
many-body matrix elements from this data. The values of M have been obtained, and these 
provide a new means to study the parity violating interaction in nuclei based on statistical 
theory 0. 

The particular formulation on which we base our paper is that of statistical spec- 
troscopy (4j. Statistical spectroscopy has proven very successful for calculations of strength 
functions of positive parity operators, but it has not been developed as fully for negative 
parity operators, especially for the case of parity violation. In this paper we will develop 
the theory in a study of parity violation for the targets measured in Ref. ||, namely 238 U 
for A ~ 230 and I04,i05,i06,i08p d for A ^ 100 . The statistical densities that we will need are 
calculated using realistic single-particle energies and a residual strong interaction adjusted 
to fit the experimental density of states in the two elements of interest here. The determi- 
nation of the smoothed strength and level densities is well established, and the theory has 
been shown to be well justified through numerical and theoretical evaluations in the early 
work [|J. The status of statistical spectroscopy in the other nuclei measured in Ref. || will 
be explored more fully in a later paper || . On the other hand, the application of statistical 
spectroscopy to parity violation raises new issues, and it is the main purpose of this paper to 
identify the new issues and show how statistical spectroscopy may be implemented in view 
of them. 

We discuss the statistical observables in Sect. II, emphasizing their advantage for study 
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of the weak PNC interaction in compound nuclei. In Sect. Ill and the Appendix, we 
emphasize the microscopic foundation of the statistical theory and also collect ingredients 
that will be used in this paper for specific calculations. The limited aim of our calculations 
of M, presented in Sect. IV, is to benchmark the theory, i.e. provide reference points for 
future calculations based on standard ir plus p meson exchange interactions with strengths 
taken from the existing literature (specifically we use the "DDH best" set, see Ref. HI). 
Since we are not adjusting parameters to fit data, we leave open the question of whether 
the data prefers some alternative parameter set f|. 

In Sect. V we discuss sources of uncertainty in such a calculation and indicate additional 
studies that would establish a more refined theoretical error analysis. Finally, in Sect. VI 
we present a summary of our paper and draw our conclusions. 



Experiments are often able to determine individual matrix elements < \l//|0|\I/i<- > of 
some operator O between stationary states > and > of the nucleus. Microscopic, 
nonstatistical calculations of such matrix elements are often quite successful. However, such 
calculations of individual matrix elements are not guaranteed to be stable. For example, 
a basic problem arises in calculating weak interaction operators O (e.g., the axial charge 
operator or the hadronic PNC interaction) within the shell-model framework. This is because 
the matrix elements, 



i.k 

which we have expressed in terms of the expansion coefficients a] of the nuclear eigenstates 
> in a many-body basis \4>i >, 



often involve strong cancellations so the predictions are sensitive to small changes in the 
Hamiltonian of the strong interaction. 



II. STATISTICAL AND NONSTATISTICAL LIMITS 




(1) 




(2) 
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It is instructive to see how this sensitive interplay comes about. In order to avoid 
the problem of diagonalizing the complete many-body Hamiltonian, one assumes that the 
important pieces of the nuclear wave functions are the low-lying components, which one 
isolates by making a truncation into a "model space" . Because the components in the model 
space correspond to just a piece of the wave function, one must make a correction, either by 
applying a renormalization to the hadronic weak interaction to account for the excitations 
out of the model space not accounted for by the effective interaction, or by enlarging the 
model space. In the shell-model theory of Ref. ||, this was done by enlarging the space. 
Enlarging the space was accompanied by corresponding changes in the coefficients a\ and 
a k K) leading to dramatic changes in the matrix elements, illustrated very clearly in Ref. ||. 
It is the conclusion of this and other investigations that because of the sensitive interplay 
among various components, the nuclear wave function must be calculated very precisely, 
especially by including multiple particle-hole (p-h) admixtures and carefully eliminating 
spurious center-of-mass motion. To avoid this strong model dependence, it was pointed out 
that the PNC matrix element could be calculated from the measured analog first forbidden 
decay of 18 Ne. 

With more computer power the individual matrix elements may of course be calculated 
in light nuclei with sufficient stability. However for heavy nuclei, where one has many more 
degrees of freedom, chaotic behavior may occur, and then the problem becomes much more 
serious. In this case eigenstates of the Hamiltonian have so many components, with ampli- 
tudes roughly equal in size and with phases that occur in a nearly random fashion, that the 
individual matrix elements become impossible to interpret in principle. Even small changes 
in the Hamiltonian (or approximation scheme) can dramatically affect the wave function and 
hence the expectation value of an observable 0. (The effective dimensionality, or number of 
significant Slater determinant components in the eigenstates, is far more stable.) Although 
a large amount of information is irretrievably lost in the transition from simple to complex 
systems in this fashion, significant physically meaningful information does remain. Indeed, 
it was the original motivation in developing the theory of statistical spectroscopy to identify 
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this information and develop methods to extract it. The chaotic behavior of the nuclear 
wave function suggests application of the central limit theorem for this purpose and leads to 
the notion that matrix elements of operators may be thought of as elements of a statistical 
distribution. 

For chaotic nuclei, one is thus led to identify as the observable the square of the matrix 
element of the PNC interaction averaged over Ni and Nk states of appropriate symmetry 
occurring in a small window of energy, 

M 2 = -±—Y,\<^i\0\^k>\\ (3) 

which can be cast into a form very similar to Eq. (|I|) using the statistical properties of the 
wave function. To see this, first use Eq. @ to rewrite Eq. (|j) as follows, 

m 2 = E E(^E«r«/) < K\o\k >< K\o\^ >* (|j:«f4)- (4) 

Since the phases of the components are random when the statistical approach is valid, we 
may write 

^- E 4 2 *4 = ^)^E I4 2 l 2 , (5) 

and if the a K (1 < k < n, where n is the effective dimensionality) in the averaging interval 
are additionally assumed to be elements of an ensemble of Gaussian distributed random 
variables, we may write for each k 2 

^EI4 2 I 2 4- (6) 

Typically, Nj and Nk are rather small in comparison to the effective dimensionality n. 
Equation (^j) is the stronger assumption. In Ref. ||, the space is partitioned with the 
different partitions being described by separate ensembles. If we adopt the weaker statistical 
assumption of Eq. (||), Eq. ([|) becomes 

m 2 = t^EEK* < ^\°\^ > 4l 2 (7) 



Comparing Eq. (0) to Eq. (|l|), we see that now the different terms in the sum over i, k occur 
with the same sign, so that the problem with delicate cancellations that would prevent 
a meaningful calculation of an individual matrix element is avoided. Thus, observables 
amenable to a statistical analysis are robust against approximations that do not change the 
overall scale characterizing the distribution of the a}. 

In actual situations where the statistical theory is applicable, only the low-lying compo- 
nents a k K of the many-body wave function are statistical in the sense of Eq. @. We look 
at a more systematic method for separating the statistical and nonstatistical behavior in 
Sect. Ill below, and find that when a more careful separation is made, we arrive at the same 
conclusion. 

III. FORMULATION OF STATISTICAL THEORY FOR PARITY VIOLATION 

The specific problem we will be addressing in trying to understand the TRIPLE data of 
Ref. [[J is how to calculate M when the operator O in Eq. (^) is the interaction V PNC , cor- 
responding to the hadronic parity- violating interaction between two particles in free-space; 
we base our numerical results in this paper on a meson-exchange form for the weak PNC in- 
teraction with the coupling strengths presented by Desplanques, Donoghue, and Holstein || 
(DDH). Although these coupling strengths were estimated by DDH from empirical non- 
leptonic hyperon decay data, they also incorporated some estimates of various contributions 
such as the strong interaction enhancements to the symmetry breaking; as a result, the cou- 
plings still contain a considerable amount of uncertainty. Because of these uncertainites and 
the fact that the DDH interaction has not been fitted to nuclear PNC data, one does not 
expect to be able to explain the TRIPLE measurements without some further refinement of 
the coupling parameters. 

In our formulation, we will be using the fact that statistical methods [ffl provide a rela- 
tively straightforward and model-independent means to average over the important compo- 
nents (i,k) in Eq. (|7]), namely those that mix chaotically among themselves by the action 
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of the strong interaction. These components correspond to the independent-particle config- 
urations formed by placing valence particles in single-particle orbitals that lie within a few 
MeV of each other (within the spreading width) outside a closed shell, in other words, to 
components within a model space. The central limit theorem is used to express the strength 
functions in terms of bivariate Gaussian functions characterized by a few parameters as de- 
scribed in Sect. III. A below. Statistical mixing of single-particle states lying farther than a 
typical spreading width (a few MeV) away from these is weak [TIJ and naturally suppressed 
in the statistical formulation. 

However, statistical mixing is unfortunately only part of the problem of understanding 
the TRIPLE data. Nonstatistical mixing of states lying outside the model space is expected 
to occur (see Sect. III.B.3), and the real question at issue in applying statistical spectroscopy 
to parity violation is how to deal in a satisfactory way with important components of the 
wave function that lie just lhu outside the model space, for example the giant resonances. 
One solution to this problem is provided by effective interaction theory |[I . 



Effective interaction theory provides book-keeping techniques for building the contri- 
bution of the (i, k) components lying outside the model space in Eq. (|^) into an effective 
operator O acting entirely within the model space. More generally, it prescribes not only 
how one can calculate microscopically various effective operators O corresponding to actual 
observables O, but also the effective Hamiltonian H whose eigenfunctions |\E^ >, 

H\^ K >=E K \^ K >, (8) 

correspond to the true stationary states \^/k >• Once the model space has been set up, 
matrix elements of the Hermitian, energy-independent operators H and O may be calculated 
diagrammatically according to the rules of degenerate many-body perturbation theory |12| . 
The connection to the observable given in Eq. ([!]) is 

< */|0|*k >=< ^/|0|^ > . (9) 

An application of effective interaction theory for the calculation of PV observables in the 
sense outlined here was followed in Ref. |13[ . 
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Once one has obtained the effective operator V FNC corresponding to V PNC by making 
use of this theory, one may apply the powerful methods of statistical spectroscopy directly to 
V r1 ^. The expression for M 2 in Eq. (0) may be expressed in the model space by expanding 
> as in Eq. (H) to obtain 

\V K >=Y / a k K \(p k >, (10) 
k 

where the sum now runs over independent-particle configurations k lying purely in the model 
space. (Note that the coefficients a k K do not necessarily have a simple connection to a k K .) 
By substituting Eqs. @ and (|H]) into Eq. (|l|), we find 

< */|0|^ >=< ^i\0\^ K >= £a? < UO\<t>k > a k K . (11) 

i,k 

Similarly, the model space counterpart to Eq. (^) is found by substituting Eqs. (|TT|) into 
Eq. @ and following the analogous reasoning leading from Eqs. (f|) to Eq. ([?]), 

M 2 = ttVEE I*? < <t>i\0\<!>k > 4l 2 - (12) 
The connection between M in Eq. ( |12|) measured in a region of excitation energy E, the ef- 



fective parity-violating interaction, and the strength-function defined below in Sect. III. A 
(with O = V PNC ) is 

M 2 = M 2 (E l7 E 2 , r 1; T 2 ) =D(E h T 1 )D(E 2 , T 2 )3 U {E 1 , T u E 2 , T 2 ), (13) 

where D(E, T) is the mean spacing between neutron resonances and T is the set of conserved 
quantum numbers labeling the states of interest (i.e., angular momentum, parity, isospin, 
etc.). 



As we have emphasized, in Eqs. (|Tl|) and ( p!2|) the sums run only over independent- 
particle configurations i, k lying within the model space. By comparing these equations we 
arrive at a conclusion similar to that found at the end of Sect. II, namely that the delicate 
cancellations that occur among amplitudes a*/ < 4>i\0\(f)k > a k K in calculations of individual 
matrix elements are absent when one calculates the mean-square matrix elements in effective 



interaction theory. However, we also see that some residual interferences may creep back 
into the statistical problem through the calculation of O, whose various terms may have 
fluctuating signs. Since these signs do not arise from statistical behavior, one may expect 
that this method of calculating mean-square matrix elements is the more robust procedure 
when statistical behavior is indicated. 

Note how the statistical and nonstatistical aspects of the many-body physics cooperate 
to produce physical observables in this formulation: the statistical behavior is confined 
to the model space and handled within statistical spectroscopy, whereas the nonstatistical 
behavior (including any l?ku collectivity that may arise from giant resonances, for example) 
is concentrated in the effective operator terms O — O. These two aspects are interwoven at 
the end when M is evaluated in Eq. fljjj) or Eq. fllcf) . 



A. Statistical Spectroscopy of French and Collaborators 

From a statistical point of view, spectroscopy in many-body systems reduces to two 
largely independent physical features, fluctuations, and second, the secular behavior of sys- 
tem dependent properties such as level densities, expectation values, and strength densities. 
Random matrix theories, first introduced by Wigner [14|, and its growing number of general- 



izations especially with regards to mesoscopic, many-body electron systems (HJ, address the 

physics of fluctuation properties. This is the most widely recognized aspect of a statistical 

approach, but is not our concern here. Instead, we shall only describe an implementation of 

statistical spectroscopy that describes system-dependent large-scale secular behaviors that 

ptvc 

arise in the construction of the strength function of V in chaotic regimes of heavy nuclei. 



Just as Weyl Jl6| was able to show that to leading order in wave number the only surviv- 
ing information in the count of modes in a cavity was its volume, generalized central limit 
theorems (CLTs) and moment methods, as developed by French and co-workers over the 
past thirty years ||, identify the analogous surviving information in a microscopic approach 
to many-body spectroscopy and provide the means to calculate secular behaviors. In this 
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context, convergence of the CLTs improves with the dilute limit of many valence particles 
occupying a much larger number of available single particle states. It should be noted also 
that experience has shown good convergence to the CLT forms even in ds shell model exam- 
ples [|j] that involve spaces much smaller and less dilute than those appropriate for medium 
and heavy nuclei. 

The matrix element M defined in Eq. (§) is an example of a strength density, 
Sq{E^ Ti, E 2 , r 2 ), which may be expressed as a trace. In statistical spectroscopy, we con- 
sider the trace decomposed into its sum over particle configurations, M = (mi,m 2 ,m 3 , ...) 
where each rrii is the number of particles in a subset of single particle levels and N = 
(Ni, N 2 , N 3 , ...). The grouping is not restricted to a particular scheme, but for the present 
discussion each subset is taken to be a j-orbital. 

In general, strength densities are bivariate functions of energy. Using mainly the notation 
of French et al. J|], i.e. (( )) indicates the trace and ( ) indicates the expectation value 
((( ))/N), the expression is 

So(E l7 r i; E 2 , T 2 ) = 0P r2 5(E 2 - H)OP ri 6(E 1 - H)}) (14) 

where H is the effective nuclear Hamiltonian, Pp is the projector onto the T subspace, 
and O is, for example, the effective parity violating interaction. It is useful to perform a 



partitioning decomposition of Eq. (|Tj) 



So(E u T 1 ;E 2 ,T 2 )= J2 ((O j Pr 2 5(E 2 -H)P M ,OP ri 6(Ei-H)P M )} 

M,M' 

= £ So(E 1 ,T l ;E 2 ,r 2 ;M,M') (15) 

M,M' 

Parity is taken care of trivially since the configurations respect it and one can just restrict 
the M, M' sums to the appropriate parities. For the model spaces of the medium and heavy 
nuclei considered in this paper, only the angular momentum decomposition remains to be 
taken into account. If O were an electromagnetic transition operator that coupled different 
angular momenta, we would have to modify the angular momentum decomposition given 
below. However for the case we are interested in here, O is a parity violating interaction and 
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is a J-scalar. In Eq. ([14]), Pr 2 commutes with O, H and can be translated next to Pr 1 - The 
strength expression vanishes unless r 2 = Ti = T. There are two additional simplifications. 
For the purposes of the angular momentum decomposition relevant to the parity violation 
calculations, E\ = E 2 = E. Finally, as a consequence of statistical spectroscopy, the spread- 
ing widths are nearly constant. Therefore, Sq(E, E, T) is proportional to the J-decomposed 
density of states, where the constant of proportionality C(j\E) is the conditional probability 
density of finding the angular momentum j given E. We find, 

%(-^i> Fi; E 2 , T 2 ) = Sq(Ex, E 2 ; V) 

*C(j\E)So(E,E) 

= C(j\E) So(E,E;M,M') (16) 

M,M' 

To the level of our statistical approximations, C(j\E) is well approximated by a partitioned 
version of Bethe's form [17 with an energy-dependent spin cut-off factor, crj(E, M), evalu- 



ated as 



8iraj(E, M) J *\ 8a] (E, M) 
3a|( J E;,M)^(J 2 ) M + (J 2 ( J ff- J E;M)) M ^^ 1 + (17) 



a 



M 



where 



E M = ^2m i e i (18) 



and 



2 A/ 2\m \- (N r - m r )(N s - m s - 5 rs )m t (m u - 5 tu ) v^rrm/r 12 no^ 
^ (7V r - d rt - d ru ){N s - d st - 6 SU - d rs )N t (N u - 6 tu ) ? 

t<u 

(below we define VjF stu ) with e, the single particle energy of the i th orbital including any 
contribution to the energy from the part of the two-body interaction which transforms like 
a one-body operator. It is further supposed that the single particle basis diagonalizes the 
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full one-body part of the Hamiltonian. The V^ stu are the matrix elements of the remaining 
two-body interaction. The counting factor, [r], reduces to 2j + 1 above. 

In addition, only those terms survive in Eq. (|T6| ) in which the fundamental matrix ele- 
ments of O appear as absolute squares. We thus have 

So(E 1 ,E 2 ;M,M')=T[( Nl ) ( Nj ) I (M'|0|M> | 2 p(^i, E 2 ; M, M') (20) 

i,j ^ m i' / \ m j / 

to evaluate. If O is a two-body operator, the matrix element can be deduced from Eq. ([19]) 
with the {r, s, t, u] sums suppressed and if not, from the appropriate analogous equations for 
other rank operators. The theoretical expression for p(E±, E 2 ; M, M'), which can be derived 
assuming statistical properties of many-body wave functions and has been confirmed in 
numerical tests, is that of a unit normalized bivariate Gaussian, 

p(E u E 2 ; M, M') = 1 7==^ exp (- — 1 { 

2£{Ei — E M )(E 2 — E M >) (E 2 — E M /) 2 ] 

C"MO"M' 0"m' J 

which to specify completely requires the evaluation of two centroids, -Em, -^m', two variances 
cr M ,cr M , and a normalized correlation coefficient, £. In principle, these moments are given 
by 

(0Pm>OP m )) Z = ((6 j HPM>OHPM))/(S 00 a M a M ,) 

((O j P M ,OHP M ))/S 00 E M , = ((O j HP M ,OP M ))/S 00 
(0P M ,OH 2 P M ))/S OO a m , = (0H 2 P M ,OP M ))/S OO 

(22) 

However, the centroids and variances can be replaced by the simpler expressions given in 
Eqs. ([l8|Jl9D with almost no loss in accuracy of the approximation. 

The correlation coefficient £ is quite important. As it approaches zero, strength is dis- 
tributed uniformly in the sense that the many-body mean square matrix element is a constant 
independent of {E\, E 2 ). In the opposite limit of £ — > 1, the mean square matrix elements of 
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5*00 — 
Em = 

M = 



O are concentrated along the diagonal and their behavior approaches being inversely pro- 
portional to the level density; see ahead. Although, £ depends on (M, M'), the fluctuations 
from configuration to configuration are small and the errors made in using a single averaged 
value negligible compared to other approximations. In fact, we use Gaussian ensembles with 
operators of fixed particle rank to evaluate £ since the strength density is not extremely sen- 
sitive to small variations. In the dilute limit of N (number of single-particle states) — > oo, 
H a /c-body operator, O a fc'-body operator, and m valence particles, 



The finite N expression is more complicated but analytically known for a Gaussian ensemble. 
It is this expression that we use and which can be found in Ref. ||. It is not sufficiently 
illuminating to give the expression here. Typical values for £ in medium and heavy nuclei 
are ~ 0.7-0.9. 

The criterion for choosing an appropriate model space is more complicated in the strength 
forms than for the level density since it depends on the joint contribution of (M, M'), and 
the matrix elements connecting configurations. Nevertheless, one uses a similar logic and 
the most important configurations are the same as in the level density calculation, i.e. the 
single-particle orbitals strongly mixed by the residual two-body interaction. If O connects 
to important configurations outside the model space, statistical methods must be extended 
to incorporate perturbative effects as discussed in the following subsection. 



Having made these general remarks, we next discuss the choice of model space and present 
the contributions to the effective one-body and two-body parity-violating interactions that 
we will use in this work. In the theory, there exist also three- and higher-body contributions 
to the effective operator O, but there have been no calculations of these terms so we will 
omit them. 




(23) 



B. Model Space and Effective Interactions 
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1. Model Space and Strong Interaction 



There are two practical issues with regard to the choice of the strong interaction in 
statistical spectroscopy. One is the specification of the one-body part of the interaction, 
which specifies the model space and single particle energies. The other is the choice of the 
spreading interaction that is crucial for getting the level density correct. Once these two 
elements are specified, then the partitioning of the model space and the calculation of the 
bivariate densities in Sect. III. A are completely specified. 

By "spreading interaction", one means that part of H that spreads the configurations 
across the eigenstates. To determine it from a given H requires removing those parts respon- 
sible for shifting the centroids of levels. In the language of Ref. [|J , the spreading interaction 
is the irreducible rank-2 part of H. The shell-model or mean-field contribution, which deter- 
mines the single-particle spectrum, is the rank-one contribution of H. Note that a somewhat 
different definition of one- and two-body parts of effective interaction is used in microscopic 
shell-model studies, where the single-particle energies are often taken to be those appropri- 
ate to the beginning or end of the shell. The residual interaction for multiple numbers of 
particles in a shell can then both spread levels and shift their centroids. One must keep 
this distinction in mind when comparing the shell-model and statistical approaches to avoid 
confusion. 

In our applications of statistical spectroscopy, we choose both the mean-field and spread- 
ing interactions phenomenologically. By so doing, we thus depart from a purely theoretical 
implementation of the theory. The mean-field part is identified with the shell-model poten- 
tial. By taking the single-particle energies to be those corresponding to such a phenomeno- 
logical potential, a smooth A dependence is imposed. The spreading interaction is adjusted 
phenomenologically to reproduce the level density in the excitation region of nuclei of inter- 
est where the weak spreading width is to be calculated. The energy dependence of the level 
densities for nuclei measured by TRIPLE have been extensively studied and experimentally 



tabulated [13], and improved level densities have been obtained also in the region of interest 



15 



for M in Ref. ||. We used these results for our numerical studies. 

The choice of model space is determined mainly by the considerations discussed in 
Sect. Ill, namely that the chaotic mixing occur within the model space. Since the strong 
spreading width of states is on the order of 2 MeV, a model space spanned by single particle 
levels lying within 2 or 3 MeV of each other is expected to be sufficient for those nuclei 
in which there are at least two or three particles (or holes) in the valence space for both 
neutrons and protons. 

For the nuclei of interest in this paper, it is a relatively straightforward matter to choose 
the model space subject to these criteria by looking at shell-model level schemes. In spherical 
nuclei, the single-particle orbitals will correspond to eigenstates in a spherical potential, and 
in deformed nuclei to Nilsson levels. Since statistical spectroscopy has always been applied 
(even in regions of deformation, for example the A ~ 230 region) exclusively in the spherical 
basis, so we will do our present calculations in this basis as well, in accord with our goal to 
benchmark the theory. We estimate the importance of deformation by including it through 
an effective one-body PNC interaction (see Sect. V.B). 

In our work, the single-particle energies of the active orbits for 239 U, given in Table I 
were taken from Ref. M and for Pd were calculated pi using a Woods-Saxon potential, 



U = Vf(r) + V ls £- S ^-^-f(r) 
r ar 

f(r) = {l + expC—)]- 1 , (24) 



with realistic parameters pO 



R = r A 1/3 , r = 1.27fm, a = 0.67 fm, 
N — Z 

V = (-51 + 33r 2 — — ) MeV, V ls = -0.44V. (25) 
J\ 

The resulting energies for the daughter 107 Pd is given in Table II. The model space is spanned 
by the set of all such configurations (N in number) that can be formed by putting m valence 
particles in a few (active) single-particle orbits lying in the immediate neighborhood of the 
nearest closed neutron and proton subshells (we refer to this as a Ohu space). Note that 
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there are both positive and negative parity levels included. The appearance of the intruder 
states of opposite parity is a feature particular to heavy nuclei: for the light nuclei for which 
parity violation has been considered, there are no such intruder states. 

In the following two subsections, Sect. III.B.2 and III.B.3, we discuss three different 

PNG 

pieces of V , a standard one-body piece, a standard two-body piece, and a higher-order 
correction. The latter will be seen to be peculiar to heavy nuclei, where the relative impor- 
tance of the standard one-body and two-body PNC interactions is different from what it is 
in light nuclei. 



2. Lowest- Order Effective One- and Two-Body PNC Potentials 

The standard two-body piece of the effective PNC interaction in the nucleus Vg^° 
that we use is related to the free-space two-body PNC interaction as discussed in Ref. |]J and 
in the Appendix. The standard one-body interaction Vg t ^ C ^ is, in physical terms, that part 
of the effective PNC interaction acting on nucleons outside the core (the valence nucleons) 
that originates from the core nucleons. Since Vg^ C is just an average of Vg t ^ C ^ over the 
core of the nucleus 0, we have a relatively simple connection between the two. Referring 
to the Fermi gas model (see Eq. (20) of Ref. JT[|), we obtain 

Vsu^^-EiCo + C^a.-pu (26) 
Po i 

where p is the density of the core and po is the central density in nuclei. The coefficients Co 
and C\ are given in Eq. (|39]) of the Appendix. For V^ C ^ as well as Vg t ^ C ^ we retain only 
the contributions from F n and F , the weak 7T-NN and (isoscalar) p-NN coupling constants, 
respectively, since these are the dominant contributions from the DDH parameter set. 

Next, we want to discuss the higher-order pieces of this interaction. The situation can 
become quite involved, since higher-order pieces of the effective one-body and two-body 
parity violating interactions can be formed from both Vg t ^ C ^ and V^ C ^ . It is our aim 
in this paper to limit attention to issues that have been discussed already in the literature, 
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so we will not enumerate the other terms in any systematic fashion. This limits us to only 
one type of higher-order term, discussed in the next subsection. 



3. Higher- Order Contributions Arising from Vg t ^ C ^ 

One of the differences between parity violation in light and heavy nuclei is the relative 
importance of single-particle transitions in PNC observables. Large shell-model calculations 
(see for example Ref. pl|j2^| ) indicate that the PNC transitions in light nuclei are largely 
single-particle in nature |l],|23|]. However in heavy nuclei, because of the appearance of 



intruder states of opposite parity in the Ohto model space, V can connect states within 

this space, and many-body mixing plays a much more important role. The fact that there 
are generally no active single-particle states that are coupled by Vg t ^ C ^ (i.e., that there 
are no opposite-parity single particle levels with the same total angular momentum inside 
the Ohu model space in a heavy spherical nucleus, even though there are intruder states) 
means that V std acts indirectly. Its contribution to the effective PNC interaction may 
be enumerated in second- and higher-order perturbation theory; in the theory of this paper, 
these are the major contributors to the difference of the operators — defined in Eq. @ 
and accounts fully for the effect of Vg td C ^ on M. In a deformed nucleus the situation is 
somewhat different, see Sect. V. 



The leading perturbative contribution of Vg td C ^ to V PNC has been discussed in Ref. p4 



In an effort to use physical ideas to identify the important correction, Auerbach [25] devel- 



oped the idea of "doorway" states for the PNC interaction. In this approach, the mixing of 
Vg t ^ C ^ is mediated by the isoscalar and isovector 0~ spin-flip giant resonances. The dis- 
advantage of Auerbach's formulation is that it is highly phenomenological. A microscopic 
expression V™^ 2 ^ for the doorway contribution to the effective two-body PNC interaction 
was obtained in Ref. [|13|] by identifying the doorway state as a collective one-particle one- 
hole 0~ phonon in the Tamm-Dancoff approximation (TDA). The result brought together 
the results of Refs. P3,p3, showing that the effect of the doorway is to renormalize the 
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perturbative contribution of Ref. [24j] by a factor (ujq/uj), where lvq is the energy of the un- 



perturbed 0~ state, and oj is the energy of the TDA phonon. The result for V™y consists 
of one isoscalar contribution, given in Eq. (4.9) of Ref. and two isovector contributions, 
given in Eqs. (4.10) and (4.11) of that paper. The renormalization amounts to a suppression 
of the isovector 0~ part V™ y C ^ by about a factor of about 3 with little change in the 
isoscalar 0~ part (which was found to be much smaller for other reasons). 

Desplanques has extended the study by considering the 0~ resonances in the random- 
phase approximation (RPA) f27j (see also Ref. [|26|]). This result shows that the multi- 



particle multi-hole correlations are very important. To get a nonvanishing contribution 
from the collective excitations in RPA, it is necessary to include nonlocal terms in the 
residual strong interaction. The combined effect of these two considerations is to undo 
most of the suppression arising from the one-particle one-hole TDA correlations found in 
Ref. fL3| . Desplanques found, essentially, that the doorway result changes by the following 
substitutions, 

ujq/uj — > — , isoscalar 

1 + A p 

uoq/uo — > — , isosvector, (27) 

1 + X p 

where A„ and A' are defined in terms of the Fermi-liquid parameters as 



Ar, = —G-\ — — Hn A — H*i Ho 

v 3 i 3 u ^ 3 1 15 



K-\G'i- l ^K + \H[-^H' 2 . (28) 



and where the effective mass of the nucleon M* is set to its free-space value (the dominant 
momenta in the phonon are momenta above the Fermi momentum pp). Using the Fermi- 



liquid parameters of Ref. [pq| , we calculate that A p = —0.21 and A^ = 0.27. 

We use here the RPA treatment |27| of the 0~ resonances rather than the TDA treat- 
ment. This is implemented by using Eq. (^) (see Appendix) with the replacements in 



Eq. (|27|). The isovector interaction, which is the only significant piece to this doorway re- 
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suit, continues to be suppressed in RPA, but the renormalization is about a factor of two 
smaller than it was in Ref. [13| . 



IV. CALCULATIONS 

In this section we describe our calculation of M and compare the results to the values 
of M deduced from scattering longitudinally polarized neutrons from the target nuclei 238 U 
and 104 > 105 > 106 > 108 Pci. These nuclei are good first cases to consider, because our calculation 
is relatively straightforward for proton and neutron mid-shell occupancy, and because the 
empirical values of the weak spreading width T w are rather well determined || for these 
cases. The relationship between and M is most precisely stated if we affix a subscript 
J on M to specify the total angular momentum of the CN states over which the PNC 
interaction is averaged and introduce Dj, the theoretical level spacing for s-wave levels of 
total angular momentum J in the region of the spectrum where the measurement has been 
performed. Then, 

T w = 2n\Mj\ 2 /Dj, (29) 

where the label J has been omitted from T w because of the expectation that this quantity is 
insensitive to the details of the distribution of CN levels. The resulting values for Me xp are 
given in Table III. Also shown in the table is the spacing D of all s-wave levels, which is the 
same as Dj only for J = targets. In the case of the daughter 106 Pd there are two relevant 
intermediate states, of J = 2 and 3, and in this case the quantity Dq is the combined level 
density for the two states. To determine value of Me xp in this case required making an 
assumption about the relationship between D and the two values of Dj. 

We report the results of our calculation of M based on the theory developed in earlier 
sections of this paper in Table IV. The contributions to M arise from two sources, the 
standard two-body parity violating interaction Vg t ^ C (Sect. III. A. 2) and the doorway 
piece V™y 2 ' (Sect. III. A. 3). To evaluate the matrix element M corresponding to V™y 2 \ 
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we replace Uq/u^ — > 0.78 in Table VII in accord with Eq. (p7[). Details of our calculation 
of the matrix elements of V™^ 2 ' and Vg t ^ C ^ are given in the appendix. The values of M 
corresponding to these two pieces of the interaction separately are given in the second and 
third columns of Table IV, respectively. In the case of the daughter 106 Pd, there are two 
different sets of s-wave resonances corresponding to J = 2 and J = 3 as discussed earlier; 
the individual values of Mj differ by less than 10%, so we have given only the average value 
M in this case in Table IV. The small value of M for 106 Pd reflects smaller level spacing 
seen in Table III and the near constancy of T m in the Pd isotopes. 

The fourth column of Table IV contains the value of M corresponding to the sum of Dwy 
and Std. These results must be added in quadrature to obtain the combined value of M. 
Note that in doing this one must allow for the possibility that the two pieces can interfere. 
We can easily verify from the values given in Table VII that in practice there is only a small 
interference between the Std and Dwy contributions. We find the interference is constructive 
for A ~ 230 but destructive for A ~ 100. 

In the last column of Table IV we show the empirical value of M taken from Table III. 
We see that for A ~ 100 the experimental values of M fluctuate around the theoretical 
values. As the weak spreading widths are expected to be much more stable than M (the 
theoretical spreading widths are stable to within about ±10 % ), we compare in Table V the 
theoretical and experimental spreading widths in the two mass regions considered. Given 
that the spreading widths are quadratic in M, we see from Table VI that the M in the 
mass A ~ 230 region is about a factor of 3 smaller than the empirical value and in the 
mass A ~ 100 region about a factor of 1.7 smaller. The agreement between theory and 
experiment is remarkably close and encourages further applications with the goal of refining 
the empirical determination of the weak coupling parameters. 

Earlier calculations ||29|J2^JT3[1 in the mass region A ~ 230 share some features in common 



with the present approach. Reference described the nucleus as a chaotic Fermi gas, i.e. 
the basis states of the nucleons were taken to be plane waves. The framework of statistical 
spectroscopy was used to calculate M with the strength function of the PNC interaction 
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normalized to the strength function calculated in a finite nucleus ||. When plane waves 
are used, the selection rules are much weaker than with shell-model states, and for this 
reason both Vg t ^ C ^ and Vg t ^ C ^ contributed directly to M. In the present work V^ C ^ 
contributes directly only if the nuclear core is deformed (see Sect. V.B). Otherwise, V^ C(1) 
contributes indirectly through spin-flip excitations. As discussed in Sect. III.B.3, Vg t ^ C ^ 
contributes directly because of the presence of opposite-parity intruder states. In Ref. [^9j M 



was found to be several times larger than the empirical value (see also Ref. |H| ) , whereas in 
the present work it appears to be several times smaller. It would thus appear that the nuclear 
matter treatment overestimates the result of statistical spectroscopy in a finite nucleus for 
the contribution of our one-body PNC interaction for the daughter 239 U. 

Reference employed a Oku model space like ours, so Vg t ^ C ^ contributed indirectly, 



through what can be thought of as the perturbative piece of the doorway process. The 
average over states was defined by a Lorentzian function approximating the spreading width 
and normalized to an empirical estimate of the number of "principal" components of the 
nuclear wave function. Good agreement between measured and calculated values of M were 



reported. In Ref. pi the indirect contribution of Vo t ^ C dominated the contribution of 



Vg t ^ C< " 2 \ and the latter was neglected. In our work, the matrix elements of V^ C ^ are 
similar to those in Ref. |24| (note that a different representation of the two-body PNC 
interaction was used in our work). However, in contrast to Ref. ||24|| , the matrix elements 



of Vg t ^ C< ^ are comparable to the sizes of the matrix elements of our V™^ 2 \ as can be 
inferred from Table IV. Although we have used different parametrizations of the effective 
strong interaction, we do not believe that this can fully account for the large differences we 
find. 

It was suggested in Ref. |24j] that the contribution of Vg t ^ c ^ to M, relative to that of 



v std C{2 \ would follow an Al/3 law - We 

see from Table IV that the relative value of the 



doorway contribution does increase with A as predicted in Ref. but a more detailed 
estimate suggests that it increases by more like a power of A 1 / 2 . The existence of the different 
rate of increase of the standard and doorway contributions shows that the A-dependence 
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of the measurements may provide an empirical means to determine the size of the doorway 
contribution. 

V. DISCUSSION 

In this section we will address the implementation of statistical spectroscopy and various 
sources of uncertainty in the calculations of the previous section. Most of the uncertainties 
we discuss are amenable to further analysis in large-basis shell-model studies, which would 
be valuable to gain additional confidence in the theoretical results. 

The theory of level densities and strength distributions within statistical spectroscopy 
applies to individual systems and operators, and as such is not inherently an ensemble theory. 
The central limit theorem results and correction expansions are all based on having many 
strongly interacting particles dilutely spread among even more single-particle states. The 
inputs required for all calculations are low-order operator moments. If it were possible to 
calculate exact moments up to products of approximately, say, six operators decomposed for 
the symmetries involved, we would have an essentially "perfect" implementation of statistical 
spectroscopy. It would be extremely accurate and would work even far out in the tails of the 
distributions. Instead, we have made a number of practical compromises in order to complete 
this first implementation for the parity-violation problem. Some of the compromises can and 
should be eliminated in future, planned, improved calculations. 

Most of the compromises we have chosen are motivated by the following considerations: 
i) the necessary angular momentum decomposed moments cannot be written down analyt- 
ically; ii) even the non-decomposed moments (scalar moments) are not worked out beyond 
a product of four operators-the expressions for products of three or four operators, though 
they exist, are almost prohibitively unwieldy; and, in) the exact operators may not be 
known in the model spaces-for example, the exact, residual, strong operator in a model 
space appropriate for heavy nuclei. We have addressed these points by: i) treating the an- 
gular momentum as a statistical variable in the spirit of Bethe's ansatz for the spin cutoff 
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factor; it) use partitioning of the moments to improve accuracy instead of incorporating 
polynomial corrections to the lowest order central limit theorems-this alleviates much of the 
need for higher moments; and, in) supplement exact moments with ensemble calculations 
and nuclear modeling-the idea being that in chaotically behaving systems whole classes of 
operators deviate very little from certain ensemble results (subject to very few restrictions). 
This also reduces the need for higher-order moments. 

To be more specific about ensemble results, we have in mind two-body random operator 
ensembles which are a subclass of embedded Gaussian ensembles ||30|| . The basic ensemble 



would have random k-body matrix elements embedded in an m-particle space distributed 
over N single-particle levels; two-body {k = 2) is the most important one. Some of the 
moment simplifications that we are making in this paper, see the comments after Eq. (p2|). 
are ensemble results. Expectation value relations such as 

(OH) = (0)(H) 
{0 2 H 2 } = (0 2 }(H 2 ), (30) 

which express a factoring approximation true in systems that behave statistically, have been 
found to hold to a considerable accuracy (see Refs. ||||). Operators O or their higher 
powers, which effectively contain a piece proportional to H (or a function of H), will not 
satisfy relations of this kind. 

In the following subsections, we address issues related to the nuclear modeling we have 
done, specifically the size of the model space; deformation; correlations between O and H 
including spurious center of mass motion; and, the choice of weak coupling parameters. We 
critique each uncertainty below, in Sects. V.A to V.D. In Sect. V.E we discuss how the 
shell-model could be used to gain additional confidence in the models. All of the statistical 
spectroscopy compromises that were discussed above could (and should) be incorporated 
into the shell model testing. 
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A. Size of Model Space 



The physical motivation for choosing a OTiuj model space was discussed in Sect. III. 
However, to establish that the optimal selection of active orbits has been made, it must be 
verified that the amount of parity mixing as determined by our statistical methods are stable 
under enlargements of the model space. For this purpose, we added a spherical I/5/2 neutron 
and l(?7/2 proton orbit to the model space for 239 U. Adding these two orbitals, found that the 
calculation of M decreased by about 10% when the residual strong and weak interactions 
were chosen as described above. Note, according to the theory of Sect. Ill, the effective 
interactions should be changed as the model space is enlarged. Although we have not made 
such compensating changes, we expect that the effect is small for the small enlargements 
considered here. Our results thus indicate that the estimates based on the model space given 
in Table I and Table II are sufficiently accurate for the purposes of comparing to experiment, 
where the errors are greater than 10%. 

The reason for the stability is that the neutron and proton valence orbits are not close 
to being either fully occupied or completely empty. In the opposite situation, choosing the 
optimal model space becomes a more delicate issue. There are a number of cases measured 
in Ref. [[J in which this is true for one of the shells (for example, Nb and Cs), and in these 
cases the choice of the optimal model space requires more effort. The data show that the 
spreading width is anomalously small in these cases, and it will be an interesting test of the 
theory to be able to reproduce this trend. 

B. Deformation 

Deformation contributes in two possible ways: if the nucleus has a permanent defor- 
mation, the one-body PNC potential develops a nonspherical component arising from the 
interaction of a valence nucleon with the deformed core. Otherwise, the two-body effective 
interaction may have a contribution that arises as one of the nucleons interacts with a core 
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nucleoli, exciting it and causing a core-deformation fluctuation. A second nucleon may then 
experience a long-range interaction with the first through the resulting core polarization. 
We estimate the effect of permanent deformation in this section. 



1. Permanent Deformation 

The fact that the empirical and theoretical values of M more nearly agree for A ~ 100 
than they do in the region of A ~ 230 could point to an important contribution of permanent 
deformation, since 238 U is a highly deformed nucleus but Pd is not. It is easy to show that 
there can be non-vanishing matrix elements of the deformed piece of Vg^ C for opposite 
parity spherical orbits in which the j differ by zero or two units (and for opposite parity 
deformed orbits in which the z-projection quantum numbers Q are the same). For the 
spherical orbits in Table I there are pairs of orbits for both protons and neutrons that 
satisfy this requirement and have rather large matrix elements. Thus, for our calculations 
in the A ~ 230 region we will have direct contributions from V std even in the Ofiu model 



space, similar to the situation in nuclear-matter ]29 . 

We find an estimate for the contribution to M arising from the deformation by using 
a one-body potential similar to that given in Eq. (|2~6D and using an Harmonic oscillator 
basis, replacing p by a density having a deformed Woods-Saxon shape with radius R and 
diffuseness a. The contribution of the one-body PV potential arising from deformation is 
found by expanding the density to lowest order in the deformation parameter /?2, leading to 
the following shape, 

6p(r)/p = -^g( r )Y 20 (r) (31) 

CI 

where 

n(r) P ( r - R )/ a 

g( r ) = P{ ] —. p .. . (32) 

Taking matrix elements between the single-particle orbitals, we obtain 



26 



< n f l f 3 f m i \5V PNa V>\n i t i j< mi >= imuJ ° R ^ ( Co + C u C - d) 

a 

< n f £ f j f m i \^(g(r)a-p + a-pg(r))Y 20 (r)\n i £ i j i m i >, 



(33) 



where at p = p , Cq + C\ = 3.53- 1CT 8 and Cq — C\ = —3.15- 1CT 9 for (proton, neutron) orbits, 
and where the radial quantum number n is the number nodes in the radial wave function, I 
is the orbital angular momentum, j is the total angular momentum, m is the projection of 
the total angular momentum along the z-axis, and 

< rififjfmilhg^a ■ p + a ■ pg^Yioir^ndijim, >= (-l) mi ~^(2j f + l)(2j< + 1) 



ji 2 Jf 



'5(2£, + 1) 



Air 



Ji 2 Jf 
mi —rrij 



y/2(t f ± 1) + 1 



± 



12^ 







< n f£f I 7:(9(r)cr -p + a ■pg(r))\n i £ i ± 1 >, 



(34) 



with Jy- 



± 



2' 



Matrix elements of Eq. ( p3|) do not vanish in the model space of Table I. There are two 
non-vanishing matrix elements, one for the pair proton orbits (^9/2,^13/2) and one for the 
pair of neutron orbits (^15/2^11/2)- The matrix element of §V PNC ^ of Eq. (^) is about 
2.4 times smaller for neutrons than for protons. Averaging the square of Eq. (|33|) over the 
quantum number (noting that there are 2j < + 1 nonvanishing matrix elements), we find 
that the magnitude of the average proton matrix element is 0.011 eV. 

We have estimated M corresponding to these values by normalizing to the calculation 



in Ref. M along the lines explained in Ref. [29], which gives 



M 2 Def (keV 2 ) = 2.6a 



(35) 



where 



1.2) 2 Tr[(5V PNC Wy 



Tr[(U 2 ) 



(36) 



The factor of 1/A V arises because the trace is taken over the squares of the one- and two- 
body matrix elements within the model space; the square of the matrix elements is weighted 
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by the relative number of valence nucleon pairs, and the factor of 1.2 2 represents the relative 
normalization of the strength function for a one-and two-body operator. The two-body 
interaction U2 has no diagonal matrix elements but is otherwise modeled after the surface 
delta-function interaction pifl . We estimated in this fashion that in 238 U, M^ef = 0.004 
meV. 

Compared to the values of M arising from other sources in Table IV, we see that the 
contribution of deformation to M is completely negligible. A more quantitative calculation 
would entail a calculation using a deformed model space of Nilsson orbits, which would 
require no new theoretical development because statistical spectroscopy as formulated in 
Sect. III.B applies in either a spherical or deformed basis. 

2. Core Polarization 

Core polarization as described in the introduction to this section has been studied phe- 
nomenologically, and its contribution to the effective nucleon-nucleon interaction is known 



to act dominantly in T = 1 and be largely independent of spin [52fl. Just as in the case 



of the doorway model [13|] the effective parity-violating interaction results from the com- 
mutator with the residual strong interaction and a ■ r. This commutator does not vanish 
for the core-polarization piece of the effective nucleon-nucleon interaction, and the resulting 
contribution to the effective PNC interaction deserves quantitative study in the future. 

C. Correlations between O and H 

There are several classes of correlations that we wish to mention. One arises when O 
couples strongly to a mode of motion generated by H lying predominantly outside the model 
space. In such a case, one would miss important physics by averaging O over only those 
states within the model space. To deal with this situation properly, one must use effective 
interaction theory and consider instead of O its effective operator counterpart O. The book 
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keeping rules of the theory ensures that corrections necessary to represent the effect of this 
correlation are included in O. 

Such a correlation is quite important for parity violation, where the isoscalar and isovec- 
tor 0~ spin-flip resonances couple strongly to the one-body piece of the parity violating 
interaction, so in this case O is identified with Vg t ^ C ^\ In order to deal with this situa- 
tion in the context of parity violation, the doorway model was introduced, as discussed in 
Sect. III. A. 3. We showed that the effect of the 0~ spin-flip excitations is quite important. 

Yet another type of significant correlation occurs when the operator O and Hamiltonian 
H commute with each other, or where both commute with a third operator R. Let us 
consider the latter case, so that H and O may be characterized by the same quantum 
numbers Y that are the good quantum numbers of R. Values of the mean-square average 
of O over distinct ensembles of states, each characterized by its own value of T, would in 
general be different. Thus, it would be misguided to average O over an ensemble of states 
with mixed T, since H does not produce such mixtures. If a statistical ensemble of mixed 
T were formed, either because one was not aware of the symmetry R, or because of some 
unfortunate approximation, this ensemble would fail to properly describe the system. 

One place where such a consideration is relevant is in the description of spurious center- 
of-mass motion. Generally, the eigenstates of H may be grouped into sets distinguished 
by quantum numbers specifying the center-of-mass motion. Suppose that each such set is 
spanned by distinct individual-particle basis states. Then, any mean-square average of O 
corresponding to one set would differ from that of any other, because the corresponding 
statistical ensembles are different. For example, if one compares two ensemble averages of 
O, one in which the center-of-mass motion is pure, and one in which it is a mixture of 
several modes, one may find different results. Fortunately, this problem turns out not to be 
a serious one for us. As can be seen from Tables I and II, our Ohcu basis does not permit any 
spurious center-of-mass excitations. Thus, all matrix elements of y PNC ( 2 ) are completely 
nonspurious. As is generally the case in applying perturbation theory to include excitations 
outside the model space, some level of spuriosity is unavoidably introduced in calculating 
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yPNC{i)_ However, the nature of statistical calculations, where the observable is a mean- 
squared matrix element insensitive to relative phases, considerably reduces the magnitude 
of the problem compared to the standard shell-model case. 



D. Choice of Weak Coupling Parameters 

It is important to bear in mind that the relative sizes of M corresponding to and 
yPNC(2) ^ e p enc j U p OI1 the choice of weak coupling constants. The values of these coupling 
constants are not certain, so the relative sizes of M for the two contributions may actually 
be different (perhaps substantially so) than the values reported in Tables III-V when more 
refined choices of the coupling constants are made. As we have indicated, we use the DDH 
parameter set for our calculations given in these tables. A more detailed breakdown of the 
operators involved as well as a careful statistical study of the additional empirical information 
available from TRIPLE and from other sources should allow extraction of a consistent set 
of values for and Fq from experiment. The application of our statistical theory to the 
remaining nuclei of the TRIPLE data set requires a more elaborate study because, as we 
have remarked, in cases not examined here the valence shells contain fewer nucleons and are 
therefore more sensitive to details requiring further refinement, such as the size of the model 
space and possibly the single particle energies. 

Note that the size of the Dwy contribution in our calculations is determined predomi- 
nantly by our choice of the weak pion-nucleon coupling constant. If, instead of the DDH 
value for this coupling parameter, we would rely instead on empirical values inferred from 
18 F H, Dwy would be considerably smaller. Likewise, the value we found for the Std compo- 
nent of the force may change substantially. For example, if the measured anapole moment 
of 133 Cs is taken to determine the weak ir and p coupling constants, with F n constrained by 
the measurement in 18 F, the value of the weak p coupling constant would be substantially 



larger [ 33 1 , as would our value of the Std contribution in Table IV. The interplay between 



the role of the Std and Dwy contribution clearly depends on the region in F n - F space of 
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interest, the details of which will be reported elsewhere. 

E. Shell Model Tests 

There are a number of useful tests of the validity of the statistical model for parity viola- 
tion that could be carried out within large-basis standard shell model calculations. We give 
one example: how to model parity mixing, as manifested in M, as it develops between states 
within the Ohu model space and those lying outside it. We have assumed in this paper that 
the mixing occurs through the doorway mechanism as modified by Desplanques f27f . The 
assumption can be evaluated in a toy model (based on large but finite spaces and schematic 
interactions) as follows: On the one hand, the chosen interaction may be diagonalized nu- 
merically (thus obtaining a model-exact result) using a large-basis shell model code. This 
same problem may then be examined in a model space applying a version of effective in- 
teraction theory. Various questions may be answered by comparing results obtained from 
the exact diagonalization and from effective interaction theory. Limiting cases would be 
instructive to consider. 

For example, if excited states are initially well separated from the model space by a large 
gap, the leading terms in the perturbative expansions discussed in Sect. Ill will dominate. In 
the opposite limit of a very small gap, which corresponds to the well defined case of statistical 
mixing in the full space, the calculation again becomes simple. The more difficult question 
of theoretical interest is how to handle the intermediate cases, where the gap is neither 
small nor large. It is in this regime where the doorway model has been conjectured to be 
applicable. It would be interesting to: (1) confirm the doorway model in specific numerical 
cases. It may even be possible to motivate extensions of it valid over larger gap sizes by 
making selective summations over specific classes of diagrams to account for multiple p-h 
excitations; (2) investigate the relative importance of effective three-body PNC interactions, 
which have so far been neglected in all theoretical studies. In the most favorable outcome, 
the doorway result we have already evaluated describes this mixing over gaps comparable to 
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shell spacings. Well-posed questions in such a model could lead to a better understanding 
of the theoretical uncertainties in applying statistical spectroscopy to PNC. 



VI. SUMMARY AND CONCLUSION 

In this paper, we have made the first extension of statistical spectroscopy to the case 
of parity violation and the weak spreading width in nuclei as measured by the TRIPLE 
collaboration || . We have emphasized that statistical strength function methods are advan- 
tageous in this case because the observables are expressed in terms of matrix elements of the 
effective PNC interaction averaged over squares of wave function components in a basis of 
independent-particle model states. We have also stressed that the corresponding theoretical 
results are less sensitive to theoretical uncertainties in familiar shell-model calculations even 
though the nuclear wave functions are very complicated. 

In extending statistical spectroscopy to parity violation, we find that it is essential to 
include corrections to the underlying parity violating interaction to account for nonstatistical 
correlations. Of particular importance are the spin-flip correlations between states separated 
by nhu, with n > 1. In this regard, effective interaction theory is particularly useful, and 
the nfiuj corrections have been incorporated through an effective operator. In calculating the 
effective operator of the PNC interaction we have implemented in the RPA approximation 
suggested by Desplanques [^7j . 



We have evaluated the weak spreading width for A ~ 230 in the 238 U target and for 
A ~ 100 in the ^Aos.ioojiosp^ targets using the standard estimates for the weak coupling 
parameters of Desplanques, Donoghue, and Holstein ||. The theoretical results are in qual- 
itative agreement with the experimental results: the measured M in the Pd isotopes differs 
by about a factor of 1.7 and in U by a factor of 3. This observation supports the hope that 
ambiguities present in the values of the underlying weak meson-nucleon coupling parameters 
may be settled by using the weak spreading width measured via neutron scattering from 
compound nuclei M. 
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We have investigated the sensitivity or our results to various improvements in the theory, 
such as the sensitivity to enlargements of the model space and to permanent deformation. 
Permanent deformation seems to be a small correction, and our estimates indicate that 
results for M are stable to 10% accuracy under enlargements of the space. The reason for 
the stability is that these nuclei are in regions of the periodic table where large numbers of 
neutrons and protons occupy the shells. When, on the other hand, shells are nearly empty 
(or nearly full), choosing the optimal model space becomes a more delicate issue. We have 
indicated that the explanation of the anomalously small spreading width in the Nb and Cs 
cases [|J will constitute an interesting test of the theory. Because of the extended space 
required, these calculations require considerably more effort and have therefore not been 
considered in this paper. We have suggested additional toy-model calculations that might 
provide insight and improve confidence in the theoretical results. 
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Appendix: Matrix Elements of Effective Parity- Violating Interaction 

In this Appendix we give explicit expressions for the effective two-body PNC interac- 
tion and their matrix elements in an harmonic oscillator basis that have been used for our 
numerical calculations. 

We may represent V PNC in general as 

V PNC = ^V PNC , (37) 

a 

where V£ NC are terms depending on spin operators 6 S ^ and isospin operators 9 T{a l The 
free-space PNC interaction in the meson-exchange model is given in Table VI. Although 
there have been numerous estimates of the strengths of the weak meson-nucleon coupling 
constants in the literature, for our numerical work we use the DDH "best" values of F n and 
F Q (corresponding to the values f n = 0.454 10~ 6 and h° p = —1.14 10~ 6 ). As we stated earlier, 
we omit the other coupling constants, since they are quite small in the DDH analysis. 

The interplay between the two-body strong and weak interaction is approximately de- 
scribed [TJ by the two-body correlation function / (r) given by Ref. [J34|| . The modification in 
Vg t ^ C ^ arising from the short-range correlation function occurs by multiplying the initial 
and final nuclear wave function by 1 + f(r), where f(r) describes the suppression of the rel- 
ative wave function at short distances due to the action of the nucleon-nucleon interaction. 

PNC(2) 

Otherwise, the representation of V std is the same as that for the free-space interaction 



as given in Eq. (|37|) and Table VI. 
Our representation of V™y 2 ^ is 

V™ y c(2) = J2 eS{a) ■ v ?( r ) Vf(R)9 T(a) + J2 d S{a) ■y?(RK(r)9 T ^ a \ (38) 

«=7,9,11,13 a=15,17 

where a dependence on R, the center-of-mass position of the two nucleons, is now possible 
because of the presence of a third body (the nucleus). The definition of the various quantities 
is given in Table VII. The pieces associated with the isovector 0~ resonance are proportional 
to C\, which is dominated by F n . These, in particular Eq. (4.10) of Ref. JI3"| , are much 



more important than the isoscalar piece associated with the isoscalar resonance, which 
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is proportional to Co and dominated by Fq (the numerical values of CV and C p in Eq. (4.10) 
of Ref. ]nj are misquoted there and should rather be identified with the values of c p and 
given in Eq. (p9|), given next). These constants are given in the local density approximation 
by 

N c -Z r 



C = c s o F + c s n F w 



N C ^Z C 

Ar 



C, = c v n F n + c v p F c - c , (39) 



where c s p = 7.61 1(T 3 , c s n = 1.74 1(T 2 , c£ = 1.74 1(T 2 , and c° p = 1.3 1(T 3 , and where iV c and 
Z c are the neutron and proton numbers of the core (A c = N c + Z c ). In our evaluation of 
Vdwv i n this paper, we have taken p = p . 

Next, we give an explicit expression for the matrix elements of Vg t ^ C ^ and V™^ 2 \ We 

PJ\[Q(2) 

express the matrix element of the effective interaction V in second-quantized notation 

as 



o 



v PNC{2) = \ E (2J+ 1)* < MxM2\\vF cm \Wbfo2 > {«i^L 2 ) J (^ 5 ^) J } 

where the sum runs over all quantum numbers (/3,j), where (3 corresponds to the radial 
quantum number n, the orbital angular momentum, the spin, and the isospin of the nucleon. 
Here a^- is a creation operator, and dfs 1 j 1 is related to the corresponding annihilation 
operator by dj m = {—iy +rn aj^ m . The density matrix corresponds to initial and final pairs 
of particles (of total angular momentum J) coupled to A J = 0. The matrix element is 
antisymmetrized 

< fal\M2\\Vr° m \Wl\&32 > = < fall M2\^ NC ®\\ fail M > 

< ^ i3ih\\V P j NC{2) 11(3^/3^ >, (41) 

and reduced in angular momentum only (we calculate separately matrix elements for neutron 
pairs, proton pairs, and neutron-proton pairs). 
A general expression tor V j is 
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(42) 



where the functional dependence of as, «r, k\, and k 2 on a as well as the values of v^ x , 
V£, S^, and are given in Tables VI and VII for the standard and doorway pieces of 

the effective interaction, and we have made a Moshinsky transformation to the relative and 
center-of-mass coordinates ((n£, NL, L\n\£\, n 2 £ 2 , £) is a Moshinsky bracket f35|). This is a 
generalization of the results given in Ref. [3£|. Our angular momentum conventions follow 



Brink and Satchler [j37| . We have checked our results against various special cases, and two 
codes have been written independently to check the calculation of Vg t ^ C ^ . The Vf stu of 
Eq. fllTf ) is equal to the Yj NC ^ j n Eq. (pEOP multiplied by \pl when the quantum numbers 



of particles 1 or 2 are equal in the initial or final states. 
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TABLES 

TABLE I. Ohoj model space for mass region A ~ 230. Note that we assume neutron (N) 

closure at N = 126 and proton (P) closure at Z = 82. The single-particle energies are given for 
the daughter 239 U. 



Particle 




Orbit 






fn i) 


Parity 


P 




ho /o 

y/ z 






(0,5) 




P 




il3/2 






(0,6) 


+ 


p 

J. 




1 7/2 










N 

IN 




in/2 








i 
\ 


N 




jl5/2 






(0,7) 




N 




§9/2 






(1,4) 


4- 


N 




d 5/2 






(2,2) 


+ 


TABLE II. 


Ofruj model 


space for 


mass region A 


- 100. 


Note that 


we assume neutron (N) 


closure at A = 


50 and prot 


an (P) closure at Z = 28. 


The sin 


gle-particle 


energies are given for the 


daughter 107 Pd 














Particle 


Orbit 




M 


Parity 


Orbital Energy (MeV) 


P 

1 


§9/2 




(0,4) 


+ 




7 Q4 


P 


Pl/2 




(1,1) 






Q D4 


P 


Lj/2 




(0,3) 






-11.2 


P 


P3/2 




(1,1) 






-10.7 


N 


d 5 /2 




(1,2) 


+ 




-8.44 


N 


g7/2 




(0,4) 


+ 




-7.34 


N 


Sl/2 




(2,0) 


+ 




-6.63 


N 


d 3/2 




(1,2) 


+ 




-5.96 


N 


hll/2 




(0,5) 






-4.85 
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TABLE III. The empirical weak spreading width T^ xp [3] (preliminary) for the daughters 239 U 


and ±uj ' ±uu > ±u <' 1u »Pcl. Also shown are the neutron number IN, the proton number 


z, 


the theoretical 


level spacings D, and the value of M calculated according to Eq. (|29|). 






Nucleus (N,Z) T^ x p (10 7 eV) D (eV) 




M Exp (meV) 


239 U (146,92) l-35lo:eI 20.75 (20.9) 




fi7+ - 24 
U - D/ -0.16 


105 Pd (59,46) l-40io!99 165 ( 165± 61 ) 




o o+2.4 
z - z -0.9 


106 Pd (60,46) 0.086±g;g|| 13.1 (13.3± 0.7) 




20+ - 10 


107 Pd (61,46) 0-18ioa5 159 ( 159± 24 ) 




n 7q+°- 88 

u - ' y -0.36 


109 Pd (63,46) 0-8llo.i 159 ( 159 ± 24 ) 




1 fi +2 -° 


TABLE IV. Theoretical values of M for the effective parity- violating interaction. 


Contributions 


are shown separately for the standard ( Std) and doorway (Dwy) pieces of the two-body interaction. 


A comparison of the experimental value of M given in Table III is also shown. 






Nucleus M S td (meV) M Dwy (meV) M Std+Dwy (meV) 




M Exp (meV) 


239 U 0.116 0.177 0.218 




0.67±£* 


105 Pd 0.70 0.79 1.03 




2.201H 


106 Pd 0.304 0.357 0.44 




20+ - 10 


107 Pd 0.698 0.728 0.968 




n 7O+0.88 
u -' J -0.36 


109 Pd 0.73 0.72 0.97 




1 fi +2 -° 


TABLE V. Comparison of theoretical weak spreading width to experiment. 


The value of T^ xp 


for Pd is obtained by statistically combining the results for all values given 


in 


Table III; the 


theoretical value of T^ he is, to within a few percent, the same for all Pd isotopes. 




Nucleus T™ e (10 7 eV) 




T^(10 7 eV) 


239 U 0.143 




± - c>o -0.64 


Pd 0.40 




1 7Q+1-65 
L - 10 -0.84 
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TABLE VI. Free-space parity violating interaction of Desplanques, Donaghue, and Holstein, 
Ref. [10]. 





(a) Definition of V™ c = 


O T( a ) S(a) . 


v a for the standard PNC interaction. 




a as otT 


T ( a ) = 


rp(a T ) QS{a) _ jg (as) 
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S (2) 


A (G lU (+) - F lU W) 
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(b) Definition of operators in Table VI(a). Note that the momentum operator p is defined 
as twice the relative momentum, p = Pi — P2- 



Operator Operator 



Spin: 




Isospin: 




s<u 


-(1)- -(2) 


-(0) 


1 




-(1)+ -(2) 


HO 


r 2 (l)-r 2 (2) 


S (3) 


% <r(l)x tr{2)/y/2 


,(2) 
HO 


r 2 (l) + r 2 (2) 


Space: 




-(3) 
r 10 


i[r(l) x t(2)]JV2 


u(r) ( - } 


[P, m ] 


f (4) 
6 00 


-r(l) ■ r(2)/V3 


u(r) (+) 


{p,M 


/6) 
r 20 


[3r 2 (l)r 2 (2)- r(l) • t{2)]/>/6 


i*o (r) 
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TABLE VII. Effective interaction arising from the doorway contribution as formulated in 
Ref. [18]. Some of the notation is defined in Table VI. Note that the doorway contribution is 
defined for odd values of a only. 



a 


a s 


a? 


h 


k 2 


7 


3 


4 


l 





9 


3 


2 


l 





11 


2 


3 


l 





13 


3 


4 


l 





15 


1 


3 





1 


17 


3 


1 





1 



6 T(a) = T (a T ) QS(a) = g(a s ) y a^ ( fm -l) Vg 2 (R) 



^^(Aio + AnJdtiJ M e -M^ i 

^^(Aoi-An)^? S« ^e"" ^ 



A 00 = -53.9 MeV-fm 3 

A i = 200 MeV-fm 3 

M = 4.758 fm _1 (nucleon mass) 



An = 239 MeV-fm 3 

Aio = 59.8 MeV-fm 3 

H = 3.897 fm _1 (p meson mass) 
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